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Abstract 

Small amplitude, long waves travelling over the surface of a current-carrying metal 
thin film are studied. The equation of motion for the metal surface is determined in 
the limit of high applied currents, when surface electromigration is the predominant 
cause of adatom motion. If the surface height h is independent of the transverse 
coordinate y, the equation of motion reduces to the Korteweg-de Vries equation. 
One-dimensional solitons (i.e., those with h independent of y) are shown to be 
unstable against perturbations to their shape with small transverse wavevector. 
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1 Introduction 



The motion of the free surface of an incompressible fluid has been the subject 
of much research, from the 19th century to the present day. Of the many 
fascinating phenomena that have been investigated, soliton propagation is 
certainly among the most intriguing [1,2]. 

Solitons were first observed in 1834 by John Scott Russell on the surface 
of a narrow, water-filled canal [3]. It was not until 1895, however, that it 
was demonstrated that solitons are solutions to the Korteweg-de Vries (KdV) 
equation, the equation of motion for small amplitude, long gravity waves [4]. 

Although it is not yet widely appreciated, there are nontrivial electrical free 
boundary problems. When an electrical current passes through a piece of solid 
metal, collisions between the conduction electrons and the metal atoms at 
the surface lead to drift of these atoms. This phenomenon, which is known 
as surface electromigration (SEM), can cause a solid metal surface to move 
and deform [5-15]. The free surface of a metal film moves in response to the 
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electrical current flowing through the bulk of the film, in much the same way 
that flow in the bulk of a fluid affects the motion of its surface. However, the 
analogy is not perfect — the boundary conditions are very different in the two 
problems. 

Just as solitary waves can propagate over the surface of a shallow body of 
water, solitons can propagate over the free surface of a current-carrying metal 
thin film [9,12,13,16]. Let the applied electric field be in the x direction and 
take the z axis to be normal to the undisturbed metal surface. If the surface 
height h is independent of y, the equation of motion for small amplitude, long 
waves on the free surface is the KdV equation in the limit of high applied 
currents [16]. The "one-dimensional" solitons satisfying the KdV equation are 
protrusions that propagate in the direction of the applied electric field. Their 
velocity decreases linearly with amplitude; in contrast, the velocity of a soliton 
on the free surface of a body of water increases with amplitude. 

A soliton propagating along a narrow, shallow channel of water is stable 
against small perturbations to its shape [17,18]. The issue of whether or not 
a one-dimensional (ID) soliton is stable is more complex if it is travelling 
over the surface of a sheet of water. The equation of motion for weakly two- 
dimensional, small amplitude, long waves on a fluid sheet is the Kadomtsev- 
Petviashvili (KP) equation 



where u^. = du/dx and so forth [2]. The KP equation reduces to the KdV 
equation if u is independent of y. In the regime in which the effects of gravity 
dominate those of surface tension, cr^ = +1 and the equation of motion is 
known as the KPII equation. ID solitons are stable in this case [19]. The 
equation of motion in the regime in which surface tension is dominant is Eq. 
(1) with (T^ = —1, i.e., the KPI equation. In this case, the ID solitons are 
unstable [19]. 

The KP equation is not the only possible generalization of the KdV equation 
to two dimensions. The Zakharov-Kuznetsov (ZK) equation 



is another; it describes deviations from the average ion density in an ion- 
acoustic plasma permeated by a strong uniform magnetic field [20]. In this 
case, ID solitons are unstable against perturbations of sufficiently long trans- 
verse wavelength [21-25]. 

To what extent does electromigration-induced soliton propagation resemble 
soliton propagation over the free surface of a fluid? Consider a long, rectangu- 
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lar, current-carrying ribbon of metal. Suppose that the metal ribbon has small 
thickness, and that it is passivated along all but one edge (Fig. 1). Solitons 
travelling along the free surface of the metal film satisfy the KdV equation 
and are stable [16-18] — this is the analog of a soliton moving along a narrow 
channel of water. 



y 








Figure 1. Cross section of the metal ribbon. The passivating 
layer (shown in black) is an electrical insulator and covers all 
but one edge of the ribbon (shown in gray). Since the thickness 
d of the ribbon is small, the location of the free surface of the 
metal is to a good approximation independent of y. The electric 
field points in the x direction if the free surface is undisturbed. 



As we have seen, a ID soliton travelling over the surface of a fluid sheet may 
or may not be stable, depending on the value of the surface tension. We now 
pose the analogous question for electromigration: Is a ID soliton travelling 
over the surface of a metal thin film of infinite extent stable? Answering this 
question is the primary goal of this paper. 

To begin, we shall determine the equation of motion for small amplitude, long 
waves on the surface of a current- carrying metal film for the general case in 
which the surface height h depends on both x and y. This equation differs from 
both the KP and ZK equations, and seems not to have been studied previously. 
We then demonstrate that the one- dimensional solitons are unstable against 
perturbations with small transverse wavevector. 
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2 Equations of motion 

Consider a metal film of thickness deposited on the plane surface of an 
insulating substrate. We take the z axis to be normal to the substrate surface 
and locate the origin in this plane. A constant current flows through the film 
in the x direction, and the electric field within the metal is Eqx. 

Now suppose that the upper surface of the film is perturbed (Fig. 2). Let the 
outward-pointing unit normal to this surface be n. The height of the film's 
surface above the substrate h depends on x, y and t. The upper film surface 
will evolve in the course of time due to the effects of SEM and surface self- 
diffusion. We assume that the current flowing through the film is held fixed. 
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Figure 2. The current-carrying metal thin film. The height of 
the free surface above the substrate, /i, depends on x, y and t. 
The outward-pointing unit normal to the free surface is n, and 
the electric field far from the perturbation is Eq. 



The electrical potential $ = <l>(x, t) satisfies the Laplace equation 

V"^^ = (3) 

— > 

and is subject to the boundary condition n ■ V$ = on the upper surface and 
z ■ V$ = on the lower. More explicitly, we have 

^^{x, y, h, t) = hx^x{x, y, h, t) + hy^y{x, y, h, t) (4) 
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and 



<^,{x,y,0,t) = 0. (5) 

If the initial perturbation is localized in the x direction, we will also have 

h{x,y,t)^ho (6) 

for X —> ±00 and, furthermore, 

^^{x,y,z,t) ^ -Eo, (7) 

and 

%{x,y,z,t)^0 (8) 
for X — > ±00, arbitrary y, and < z < hg. 

We assume that the mobility of the metal atoms is negligible at the metal- 
insulator interface, so that the form of that interface remains planar for all 
time. Further, in the interest of simplicity, we assume that the applied current 
is high enough that the effects of SEM are much more important than those of 
capillarity. The surface atomic current Ja is then proportional to the electrical 

— * — * 

current at the surface. Explicitly, Ja — —uMqW^, where the u is the areal 
density of the mobile surface atoms, q is their effective charge, and M, the 
adatom mobility, will be assumed to have negligible anisotropy. When there 
is a net influx of atoms into a small surface element, the surface element will 
move. The normal velocity of the metal surface is 

Vn = qMuflVl^, (9) 



where Q is the atomic volume and V| 
[26]. Note that ht = ^Vn, where g = I 
metric tensor for the surface. Using this 
Laplacian explicitly, we have 



is the surface (Beltrami) Laplacian 
+ hi + hf, is the determinant of the 
relation and writing out the surface 



dh { 1 

^ \ —[(1 + hl)d:, - Khydy]^{x, y, h, t) 



qMuQdt {^/g 

+ dy + hl)dy - Khyd,Mx,y,h,t)'^ . (lO) 

Here is the total derivative with respect to x, and dy is similarly defined. 
(Since h depends on x, we have dx^{x, y, h, t) = ^x{x, y, h, t) + hx^z{x, y, h, t), 
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for example.) Together, Eqs. (3) - (10) completely describe the nonlinear dy- 
namics of the film surface. 



3 Asymptotic analysis 

We wish to study the propagation of a localized disturbance whose amplitude 
is small compared to ho, and whose width is much larger than Hq. To do so, we 
will use multiple scale asymptotic analysis [1] . We put h — ho + where the 
constant a is a measure of the amplitude of the disturbance and ( — ({x, y, t) 
is of order unity. Let / be the characteristic width of the initial disturbance. 
We will study the limit in which both a = a/ho and S = Hq/I are small. More 
precisely, we shall consider the limit in which a and 5 tend to zero, but a/ 6^ 
remains finite. It is in this limit that the effects of nonlinearity and dispersion 
balance, and solitons can propagate. 

To simplify the description of the problem, we will set $ = —Eqx + 4>, so that 
V0 — > for a; — > ±oo. We also introduce the dimensionless quantities x = x/l, 



y = y/l, z = z/ho, i = {\q\MEoi'^/hol)t, and z,i) = (f){x,y,z,t)/{Eol). 

Eq. (3) becomes 



which applies for < 2; < 1 + a({x,y,t) and all x and y. The boundary 
conditions are 



<t>zz + S^{<Pxx + <Pyy) = 0, 



(11) 



05 = 



for z — 0, 



(12) 



4>s = aS'^{4>x - l)Cx + OiS'^hCy for i = 1 + aC{x, y, i), 



(13) 



for X 



±oo. 



(14) 



and 



+ (pyjj for < 5 < 1 and x — > ±00. 



(15) 



Finally, the equation of motion for the free surface of the film is 




I 



I 



(16) 
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where z — 1 + aC{x, y, i) and ag = q/\q\. Note that g — 1 + a'^5'^{(~ + C|) 
We next introduce the scaled variables 



[x + aqt) , 



a 



1/2 



a 



3/2 



T = 



(17) 
(18) 
(19) 



and 



(20) 



(^, 77, r) is a moving coordinate system that translates with velocity Vq = 
—qMEouQ/ho relative to the laboratory frame. Vq is the velocity of surface 
waves in the limit of vanishing amplitude and wavevector [9] . 

For convenience, we drop the tilde on z. Laplace's equation becomes 

+ a{ip^^ + ipr,v) = 0, (21) 

and this applies for < z < 1 + a({^,'r),T) and all ^ and 77. In terms of the 
scaled variables, Eqs. (12) - (16) are 

^z = ior z = 0, (22) 

= a'^iiiP^ - 1)Q + CM for z = 1 + aCiC, V, r), (23) 

C ^ for ^ ^ ±00, (24) 

^p^x + il^rjy — > for < 2; < 1 and ^ — > ±00, (25) 



and 



V9 



+d, 



(26) 
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for 2; = 1 + a({$^, r], r). In terms of the scaled variables, g — 1 + a^(C| + C^)- 
One advantage of introducing the scaled variables is now manifest: S does not 
appear explicitly in Eqs. (21) - (26). 

We shall now begin our analysis of the small a limit. We assume that for small 
a and fixed ^ and r, there is a solution with 

oo 

C=E"Xn(e,^,T) (27) 

n=0 

and 

oo 

V' = E«"V'n(e,^,r,z), (28) 

n=0 

where the C„'s and ipnS are independent of a. 

We next insert the expansions (27) and (28) into Eqs. (21) - (26) and equate 
terms of the same order in a. The goal will be to find a closed integral equation 
for (o. We will need to consider terms up to order a^. To prepare for this task, 
we note that to third order in a, Eq. (23) is 

1 1 

2 6 

= a^ip^ + aCi^^z - 1)C^ + aXv'^v + "\CA- (29) 
This result holds for 2; = 1. To second order in a, Eq. (26) reduces to 

V: r) + a^a^it 77, r) = (a| + 77, 1 + «(, r). (30) 

We begin by working to zeroth order in a. To this order, Eq. (21) becomes 
i^Ozz = for < 2; < 1. This implies that ipQz docs not depend on z. Applying 
Eq. (22), we see that in fact ipQz = for all z, ^, rj and r. We conclude that 
il^Q depends only on ^, r] and r, and we write 

V'o = ^o(e,?7,T). (31) 

The boundary condition at the free surface of the film gives no new information 
to zeroth order in a, as is readily seen from Eq. (29). 

We now insert the expansions (27) and (28) into Eq. (30) and simplify using 
the fact that ^002 = 0. This gives 
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+ CoV'iTjj?^ + Koni^inz + CowV'iz) (32) 
for 2; = 1. Eq. (32) is valid to second order in a. 

Armed with this result, we return to the zeroth order analysis. To that order, 
Eq. (32) reduces to 

^0^^ + Qorm — 0. (33) 

Now both 6*0^ and 6*0,, must vanish for \^\ — > 00. The only solution to Eq. (33) 
that satisfies these boundary conditions is the trivial one: ^0 is a constant for 
all ^ and 77. Because the zero of the electrical potential is arbitrary, we may 
set this constant to zero. 

To first order in a, we find that tjji is independent of z. We shall write 

A-Oiii.V.r). (34) 

Eq. (29) again yields no new information. The equation of motion (32) is to 
first order 

+ Q\r)ri — Co€ ■ (35) 

Our next task will be to write out the equations of motion to order a^. Eqs. 
(21) and (35) show that 

'4^2zz = -(^155 + ()im) = -Coj (36) 

for < z < 1. Integrating this with respect to z and using the boundary con- 
dition at the metal-insulator interface (22), we have ip2z = —doiZ. Integrating 
once again, we obtain 

V'2 = ^2 - ^Co^^', (37) 

where 62 depends only on ^, r] and r. Once more, Eq. (29) tells us nothing 
new. Eq. (32) simplifies considerably now that we know that ■^i^ = 0. Using 
Eq. (37) as well, we obtain 

<^gCoT + ^(Co^ee + Coirm) = ^2^^ + ^2r)?? " Cl?" (38) 
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Equating terms of order q:"^ in Eq. (21) yields t/j^zz + ip2^^ + 'ip2qTi = for 
< 2; < 1. Integrating this twice with respect to z and applying Eqs. (22) and 
(37), we find that 

i^3 = es- i(% + 92,r,)z' + ^(Cs«€ + Coir„)z\ (39) 

where < 2; < 1 and 9^ = 0^[^,r],T). We next equate terms of third order in 
Eq. (29). We obtain 

V'Sz + CoV'22z = -Cli + V'uCo^ + Cor,V'lr, (40) 

for z — 1. This can be simplified using Eqs. (34), (37), and (39). The result is 
Comparing this with Eq. (38), we have 

C^gCor + g(Co^e + Cor?r))^ + CoCo^ = "(^l^Co? + ^Ir^Co??)- (42) 

Since we have assumed that the disturbance is localized, ^o, 6*1^ and Oir, must 
all tend to zero as ^ ±00. 

Equations (35) and (42) are a closed system of equations that describe the 
dynamics of the film surface in the small amplitude, long wavelength limit. 
These equations represent a considerable simplification of the problem: in 
the original formulation, the 3D Laplace equation for $ = ^{x,y,z,t) would 
have to be solved subject to boundary conditions at a moving boundary, a 
stationary plane boundary, and at infinity. In the simplified description of the 
problem, a Poisson equation for 9i = 61(^,1], t) [i.e., Eq. (35)] must be solved 
subject only to the boundary conditions at infinity. 

If the film height is independent of 77, Eq. (35) reduces to ^1^^ = (0^. Integrating 
this with respect to ^ yields ^^i^ = + F, where F is a function of r. Both 
^^1^ and Co tend to zero as ^ becomes large, and hence F{t) = for all r. Eq. 
(42) may therefore be rewritten as follows: 

f^?Cor + ^Co«^« + 2CoCoc = 0. (43) 

Thus, when the film height is independent of 77, the equation of motion is the 
KdV equation, just as we found in Ref. [16]. 
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Before continuing, we shall simplify our notation somewhat. Let X — aq^, 
Y — T], T — t/3, u — 6(0, and ^ = GagOi. Equations (35) and (42) become 

ut + {djc + d^)ux + ^uux = -^(^xMx + ^^yUy) (44) 



and 

+ = ux. (45) 

The solution to Eq. (45) that satisfies the requirement that ^xx + ^yV tend 
to zero for X — > ±00 is 

00 00 

*(X,y) = ^ J dX' J dY'ln[{X -X')^ + {Y -Y'f]ux{X',Y').{46) 
—00 —00 

Insertion of this result into Eq. (44) yields a closed integral equation for u; 
finding such a relation was the goal of this section. However, we will find it 
more convenient to work with Eqs. (44) and (45) rather than this unwieldy 
integral equation. 

SEM does not create or destroy atoms — it simply moves them across the 
metal surface. The total mass, which is proportional to h{x, y, t)dxdy, 

is therefore conserved. This conservation law was not lost in the asymptotic 
analysis: Eqs. (44) and (45) imply that u{X, Y, T)dXdY is a constant. 

Indeed, using Eq. (45), Eq. (44) can be written 

ii2. + V-J = 0, (47) 

— * A A 

where V = dxX + dyY and 

- - 1 - 

J = Vux + -uV^. (48) 

— * 

Eq. (47) is the local statement of mass conservation and J is the current of 
the conserved quantity u. 

Equations (44) and (45) have soliton solutions given by 

u = 12a2sech^[a(X - - ^a^T)\ (49) 

and 

* = 12atanh[a(X - - 4a^r)], (50) 
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where a and Xq are constants. In the next section, we will demonstrate that 
these one-dimensional solitons are unstable against long-wavelength perturba- 
tions. 

If we set the right hand side of Eq. (44) to zero, we obtain the ZK equation. 
Thus, our problem differs from the ZK equation through the coupling to the 
field ^. The linearized versions of the ZK equation and Eq. (44) arc the same. 
Moreover, as we shall see, the ID solitons are unstable in our problem, just 
as they are for the ZK equation [21-25]. However, there are important differ- 
ences between the two problems: the growth rate of a perturbation of small 
transverse number k will be shown to increase as k'^ in our problem. In con- 
trast, this growth rate is proportional to |A;| for small k in the case of the ZK 
equation. 



4 Multiple scale perturbation analysis 

To prepare to investigate the stabihty of the sohton described by Eqs. (49) 
and (50), we shall transform to a coordinate system that moves along with 
the soliton, and rescale so that the soliton amplitiide is independent of a. 
Specifically, we set X — Xq = a~^{X + 4T), Y = a^^Y, T = a~^T, u = a'^u, 
and "if = in Eqs. (44) and (45). We shall again drop the tildes, and, in 
addition, replace X, Y and T by a;, y and t. This yields 

Ut - AUx + V'^Ux + ]^UUx = —{^xUx + ^yUy) (51) 

and 

= Ux, (52) 

where = ff^ /dx^ -\- ff^ /dy^. The unperturbed one-dimensional soliton is 
given by 

uq — 12 sech^x (53) 

and 

*o=12tanha;. (54) 

Note that 

Uqxxx - ^Uqx + UqUqx = 0; (55) 
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we will need this identity later. 

We now apply a sinusoidal perturbation with wavevector perpendicular to the 
direction of propagation: we put 

u^Uo + ef{x)e'''y+^* (56) 

and 

^ = ^0 + ep(a;)e*^2/+7t (57) 

and work to order e. Eq. (51) becomes 

^Lf - k'^f^ - ^uoxif - Px) = -if, (58) 

where L = d? / dx^ -\- Uq — A. From Eq. (52) we obtain 

Pxx - k'^P = fx- (59) 

We are interested in the effect of localized perturbations to the ID soliton. 
Therefore, we shall seek solutions to Eqs. (58) and (59) with the following 
properties: we require that 

/ ^ as a; ^ ±00 (60) 

and 

Pj; — > as x — > ±00. (61) 

We can obtain an approximate solution to Eqs. (58) - (61) if we restrict our 
attention to perturbations with small transverse wavenumber k and make a 

perturbation expansion in powers of k [27]. There is a technical difficulty, 
however. In a standard perturbation analysis, we would set 

f^fo + kfi + ef2 + ..., (62) 

p = Po + kpi + k'^p2 + . . . (63) 

and 

7 = A;7i + k"^^^ + . . . (64) 
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(The zeroth order term in the expansion (64) vanishes because 7 = for A; = 
[17,18].) We would then insert Eqs. (62) - (64) in Eqs. (58) and (59) and equate 
the coefficients of hke powers of k. Problems begin at the second order in the 
perturbation expansion: /2 cannot be made to vanish for both x — > +00 and 
X ^ —00. The difficulties grow still worse after this point, since /s and higher 
order terms diverge for large x. 

To overcome these difficulties, we will employ a multiple scale perturbation 
analysis [28] . Our method is closely related to the multiple scale analysis de- 
veloped by Allen and Rowlands for the ZK equation [25] . 

The basic problem is that even though k is small, x can be large. In our mul- 
tiple scale analysis, we develop asymptotic expansions for / and p. If we keep 
terms of increasingly higher order in these expansions, we obtain approximate 
expressions for / and p valid for larger and larger values of x. To begin, we 
set Xn = k"'x for n > 1, and consider a;, Xi, a;2, and so forth to be independent 
variables. Thus, / and p are considered to be functions of x and all of the Xj's, 
and for small k we have 

/ = fQ{x, xi, X2, ...) + kfi{x, Xi, X2, ...) + k'^f2{x, Xi,X2,...) + ..., (65) 



and 

p = Po{x,Xi,X2, ...)+ kpi{x,Xi,X2, ...) + k'^P2{x,Xi,X2, ...) + .... (66) 



Using the chain rule, we see that 



^ — ^+k'^+k'^'^+ (67) 
dx dx dxi dx2 



We insert Eqs. (65) - (67) and (64) into Eqs. (58) and (59) and equate coeffi- 
cients of hke powers of k. To completely determine the fiS and p^'s, additional 
constraints must be imposed. We shall require that fm/fm-i remain bounded 
as a; cxo for all m > 1 [28]. This ensures that A;"'/^ is small compared to 

Since / and dxP tend to zero as x — > oo, we demand that 

fm^^ as X —>■ OO (68) 



and 



as X ^ oo 



(69) 
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for all m > 0. Following Allen and Rowlands [25], we shall also insist that 

e^'^/m ^0 as a; ^ -oo (70) 

for all m > 0. The resulting expressions for / = /o + kfi + k'^f2 + ■ ■ ■ and 
dxP — dxPo + kdxpi + k'^dxp2 + • • ■ will tend to zero as a; — > — oo, in accord 
with the boundary conditions (60) and (61). 

If we retain n terms in each of the asymptotic expansions (65) and (66), we 
obtain expressions for / and p valid for x as large as /c^" [28]. To show that the 
ID solitons are unstable, we will need to work to order n = 3. For convenience, 
we will write fo^ = dfo/dx, foi = dfo/dxi, fo^i = d^fo/dxdxi and so on. 

Before embarking on the multiple scale perturbation analysis in earnest, we 
note that for large the problem simplifies considerably. Eq. (58) reduces 
to 

/xxx-(4 + A;2)/, + 7/ = (71) 



in this limit. Since 7 = 0{k), the general solution to Eq. (71) is a linear 
combination of e^^^, e~^'^^ and e^^^, where 

pi = 2-^7 + 0(^:2), (72) 

P2 = 2 + ^7 + 0(^2), (73) 

and 

p,^^^ + 0(k'). (74) 

Recall that / must vanish for a; ±00. Therefore, for large positive x, 

f ^ A+e-P^"" + B+eP^^, (75) 



while, for large negative x, 

f ^ A_e^^^ + 5_e^'3^. (76) 



Note that if the real part of is positive, 5+ must vanish, and, similarly, B_ 
must be zero if the real part of ps is negative. Eq. (76) will later provide a 
check of our analysis. 
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4-1 Zeroth order 



To zeroth order, Eq. (58) becomes 



where 



dxLfo - ^uoxifo - dxPo) = 0, (77) 



^^^ + ^0-4. (78) 



[Note that L and L differ: the partial derivative with respect to x appears in 
Eq. (78), whereas it is the total derivative with respect to x that appears in 
the definition of L.] Eq. (59) yields d^po — d^fo, which may be integrated to 
give dxPo = fo + Co, where Cq = Co{xi,X2, ■ ■ ■) is independent of x. Eqs. (68) 
and (69) show that Cq must vanish. Hence 

Pox = fo (79) 

and Eq. (77) reduces to 

dxLfo = 0. (80) 

This is a third order ordinary differential equation for /o; the general solution 
is 

fo — ho sech^x tanh x + AoX + Bo cosh^ x, (81) 

where 

X = ^ (3x sech^x tanh x — 3 sech^x + 1) . (82) 

ho, Ao and Bo are independent of x, but they may be functions of the x„'s 
with n > 1. Because fo must tend to zero for x — > +oo, the coefficients Ao 
and Bo must vanish. We conclude that 

fo = hoFo, (83) 

where Fo = sech^xtanhx. Eq. (79) now gives po = —^hosedi^x + Cq, where Cq 
is independent of x. Choosing the potential zero appropriately, we can arrange 
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to have Cg = 0, and then 



Po = — -/losech^x. (84) 



4.2 First order 



We next work to first order in k. Using the identity LFq = 0, Eq. (58) yields 

d^Lfi = -2/o:rxl + ^UOxifl " Plx " Poi) " 7i/o- (85) 



Eq. (59) becomes pi^x + 2poxi = fix + foi- Using Eqs. (83) and (84) and 
integrating with respect to x, we find that 

fi - plx - poi = ci, (86) 



where Ci = Ci(xi,X2, . . .) is independent of x. This allows us to simplify Eq. 
(85); we find that 

^xLf^ = -2hoiFoxx - iliho + 12ci)Fo. (87) 

Eq. (87) has the solution 

fi = -hoixFo + ^71^0 + ^ci + Ax + Bi cosh^ x, (88) 

where Ai = Ai{xi, X2, ■ ■ ■) and Bi = Bi{xi, X2, ■ ■ ■) do not depend on x.Q 
Since fi must tend to zero as x — oo, the coefficient Bi is zero. Eq. (88) is 
therefore 

3. 3. ,2l.l .1 



/i ={^Ai- hoij xFo - -Ai sedi'x + -A^ + — 71/10 + -Ci. (89) 



Now /1//0 must remain bounded as x — > 00; as a result 

Ai = ^V, (90) 



^ Strictly speaking, a term proportional to Fq should appear on the right hand side 
of this equation. However, if this term is included, the sole effect is to alter / and 
p by an irrelevant constant prefactor. Therefore, without loss of generality, we may 
omit this term. 
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ci = - ^7i/io, (91) 



and Eq. (89) reduces to 

/i = —hoi sech^x. (92) 

Using Eq. (90) to eliminate Ai from Eq. (91), we find that 

hoi + ^71^0 + ^ci = 0. (93) 

We now turn our attention to Eq. (86). Using Eqs. (84) and (92), we obtain 
pix — —^hoisech^x — ci. The boundary condition (69) is satisfied only if ci = 0; 
Eq. (93) therefore becomes 

hoi = -^71^0- (94) 



We conclude that 



I 

fi = -7i/iosech^a; (95) 

o 



and 

Pix = ^71^0 sech^x. (96) 



We integrate the latter equation with respect to x and choose the potential 
zero so that pi — > as a; — > oo. This yields 

Pi = ^7i/io(tanhx - 1). (97) 



4-3 Second Order 



Equating coefficients of terms of second order in k in Eq. (58) and using the 
identity LFq 0, we have 



dxLf2 — — '2'flxxl — 3/oa;ll — 2/oxa;2 — Lfu + fox 

+ 7,Uox{f2 - p2x - Pu - P02) - 7i/i - 72/0- (98) 
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Eq. (59) yields 



f2x - P2xx - Pixi - Pox2 = ^/iosech^x; 



(99) 



this result has been simphfied using the fact that Ci — and Eqs. (79), (84), 
and (86). Integrating Eq. (99) with respect to x, we obtain 

/2 - P2x - Pu - Po2 = ^ho tanh x + C2, (100) 



where C2 = 02(0:1, 0:2, • • •)• We now insert Eqs. (83), (84), (94), (95), (97) and 
(100) into Eq. (98). This gives 



a,L/2 = (1 - ^7i) hoFo, - 2ho2Foxx - (72/^0 + 12c2)Fo 

7^^^^°"^^° ~ e/io-Potanhx. (101) 



The last two terms on the right hand side of Eq. (101) must be rewritten 
before we can proceed further. To begin, note that Eq. (55) can be integrated 
to yield 

ul^Suo + ^SFoa,. (102) 



Using this identity, it can readily be shown that 

Luo ^ 2AFo^ + 4uo (103) 

and 

36Fo tanh a; = uo - 12Fo^. (104) 

Eq. (101) now becomes 

dxLf2 = (3 - ^^1) ^o^ox - 2ho2Foxx ~ ^ + ^^1) ^0^° 

-(72/io + 12c2)Fo. (105) 

Eq. (105) has the solution 
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/2 = - ^ + ^7?) ho cosh X sinh x + ^ (s - ^'y^^ tanh x 

+ ^(72/^0 + I2C2) - ho2xFo + Aihox + B2 cosh^ x, (106) 

where ^2 = A2{xi, X2, ■ ■ ■) and B2 = B2{xi, X2, . . .) do not depend on x.0 

Because f2 must satisfy the boundary conditions (68) and (70), the coefficients 
of cosh X sinh X and cosh^x in Eq. (106) must vanish: 

9 128 , , 

7i = (107) 



and 
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B2 = 0. (108) 



Eq. (106) now reduces to 



/2 = ^/iotanhx + ^(72/^0 + 12c2 + 6^2/^0) + (^^Aho - ho2^ xFq 
3 

--A2hosech.^x. (109) 



/2 must tend to zero as x — > cxo, and therefore 



/44 

— + 72 + 6A2 ) /lo + 12c2 = 0. (110) 



Since /2//1 must remain bounded as x ^ 00, we must also have 

h02 = ^A2ho. (Ill) 



All of the terms on the left hand side of Eq. (100) tend to zero for x — > cxo, 
and so 

C2 = -^ho. (112) 



Eliminating C2 from Eq. (110) using this result and A2 using Eq. (Ill), we 
obtain 



1 /7 1 



'-02 



4 V5 2 



+ 7^72 h. (113) 



^ Footnote 1 again applies. 
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Eq. (109) can now be simplified using Eqs. (110), (111), and (113). We find 
that 

f^^l{l + ^72) /iosech^x + ii/io(tanhx - 1). (114) 



Notice that although tends to zero as a; — > — 00, /2 itself does not vanish 
in that limit. 

We can now determine p2x using Eqs. (84), (94), (97), (112), (113), and (114) 
in Eq. (100). The result is 

P2x = ^(t + 7:^2) ho sedi^x - ^/io(tanhx - 1). (115) 



Integrating this and choosing the potential zero to be at a; = +00, we have 

-f- 

8 V5 



1/7 1 

= K + 0^2 ) /io(tanhx - 1) - -/io(lncoshx - x + In 2). (116) 



4-4 Third Order 



Using the identity LFq = once again and equating coefficients of terms of 
order in Eq. (58), we have 

dxLfz +dx^Lf2 + dx2Lfl + 2f2xxl + 3/ia:ii + 2fixx2 + 2/oxa;3 + 6/oa;12 
+/0III — fix — /oi + 12Fo(/3 — Pix — P2I — P12 — P03) 

= -(71/2 + 72/1 + 73/0)- (117) 
From Eq. (59), we obtain 

Pixx + 2p2xi + 2pix2 + Pill + 2pox-3 + 2poi2 " Pi 

= /03 + /l2 + /21 + /3x. (118) 

This result can be simphfied using Eqs. (79), (86), and (100), and the values 
of ci and C2. This gives 

fzx — P3xx — P2xl — Plx2 — P0x3 — 0, (H^) 

and, integrating, we obtain /s - p^x - P21 - P12 - P03 = C3, where C3 = 
cs{xi, X2, . . .). Now po, pi, and p2 tend to zero as x — > +00, and /a and psx 
must vanish in this limit as well. Thus, cs — and 

fs - P3x - P21 - P12 - P03 = 0. (120) 
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We now return to Eq. (117). Using Eq. (120) and our explicit expressions for 
7i, ^01, ho2, /o, /i, and /a, we obtain 



11 
20 

13 
60 

+2/lo3-foxx- 



+ 



7i/io(tanha; — 1] 
7i + 73) ^0^0 + 



:i2i) 



The solution to this equation is 



/a = ^ (^72 - — ) 71^0 cosh xsinhx 



11 , 
+ 80^^^° 



-Mo(lncosha; — x — 1) 



if \ 
+3Fr) —x^ — / In cosh x' dx' + x — In cosh x 

1 /251 



'24V60-^^-^^j+T6lr^ 
+hQ3xFo + A^hox + -B3 cosh^ x, 
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7i/io tanhx 



(122) 



where and are independent of x.[^ Applying the boundary conditions 
(68) and (70) with m = 3, we see that vanishes and 



^^=15- 



(123) 



Eq. (123) is the principal result of this section. We have now shown that the 
one-dimensional solitons are unstable. 

We must still show that our approximate solution can be given the correct 
asymptotic behavior. We require that /s — > as x — > +00 and that fs/ f2 
remain bounded in this limit. This yields 



11 



80 



71^0 



1 ^ 
--Mo ln(l + e-2^') + 3Fo J ln(l + e-^^'') dx' + ln(l + e 



-2x\ 



37 , , , If 169 \ , 

■— 7i/^o(tanhx - 1) + - (^73 + —71 j hoUo- 



(124) 



^ Footnote 1 again applies. 
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If X is large in magnitude and negative, 



A + 0(1). (125) 



vanishes for x ^ — 00, but /s does not. 

Suppose that x is large and negative, k is small, and kx is of order 1. Dropping 
terms of order e^^, we have /o//io = /i/^o = and /2//io = —11/15. Eq. (94) 
shows that 

/lo = exp {-\lix^ hi, (126) 



where hi = hi{x2,X3, . . .). In the limit we are presently considering, hi is 
approximately constant. Using these results and Eq. (125), we obtain 



/ = /o + A:/i + A;V2 + ArV3 + ... 
= -^hi^ exp (-^A;7ia;) (l + ^kjix^ +k^U + ... (127) 

Retaining terms up to third order in k, this may be written 

/ = -—hik'^ exp (]jx) . (128) 
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This is consistent with Eq. (76), and so / does indeed have the correct asymp- 
totic behavior in the x — > —00 limit. 



5 Conclusions 



In this paper, we found the equations of motion for small amplitude, long 
waves on the surface of a current-carrying metal film. In dimensionless units, 
the equations are 

1 1 

UT + (<9i + dl-)ux + -uux = --{"^xux + ^yMy) (129) 



and 



^xx + * 



YY 



Ux- 



(130) 
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These equations represent a generalization of the KdV equation to two dimen- 
sions, and seem not to have been studied previously. 

ID solitons traveling over the surface of a current-carrying metal thin film were 
shown to be unstable. The growth rate for perturbations of small transverse 
wavevector k is 

Rej = + 0{k'^). (131) 
15 

For the ZK [21-25] and KPII [19] equations, the growth rate Re 7 increases 
much faster than this as k is increased: it is proportional to \k\ for small k. 

The linear stability analysis carried out in this paper only applies shortly 
after the ID soliton is perturbed. What is the nature of the subsequent time 
evolution? For the ZK equation, it is known that if a ID soliton is perturbed, 
it will break up into a series of "cylindrical" solitons [29]. These solitons are 
propagating humps, and are stable. Numerical studies of Eqs. (129) and (130) 
are currently underway, and will reveal whether a ID soliton propagating over 
a current-carrying metal surface breaks up in a similar fashion after it has 
been perturbed. 
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